Josephson physics mediated by the Mott insulating phase 
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We investigate the static and dynamic properties of bosonic lattice systems in which condensed and 
Mott insulating phases co-exist due to the presence of a spatially- varying potential. We formulate a 
description of these inhomogeneous systems and calculate the bulk energy at and near equilibrium. 
We derive the explicit form of the Josephson coupling between disjoint superfluid regions separated 
by Mott insulating regions. We obtain detailed estimates for the experimentally-realized case of 
alternating superfluid and Mott insulating spherical shells in a radially symmetric parabolically- 
confined cold atom system. 



An important and generic situation presented by 
many-body quantum phenomena is that of competing 
states of matter co-existing in spatially-separated re- 
gions within a given system due to the presence of in- 
homogeneities. Transport properties of diverse systems 
such as the quantum Hall system, metal-insulator com- 
pounds, high T c superconductors, and more recently, 
cold atomic gases, are determined by the tunneling of 
carriers between conducting regions that are embedded 
within insulating regions^lj. Crucial to understanding 
thermodynamic and transport features of such systems 
is the manner in which conducting regions couple to one 
another through the insulating regions. Classic exam- 
ples of Josephson coupling in superconductors and cold 
atoms rely on an externally- imposed potential barrier be- 
tween condensed regions [21 E]. Here, on the other hand, 
we explore systems of bosons in which condensed (su- 
perfluid) regions exhibit Josephson physics mediated by 
Mott-insulating regions of the same bosons. This model 
should be germane to a diverse range of systems, in par- 
ticular, granular superconductors and high T c materials 
where Cooper pairs can be treated as the bosonic de- 
grees of freedom|3], and trapped cold atoms in optical 
lattices where the atoms are bosons [5]. Through an ex- 
plicit description of these phases in terms of microscopic 
parameters, we are able to go beyond phenomenological 
treatments for obtaining transport co-efficients in these 
systems pQ. 

Towards understanding this physics of co-existent 
phases, we study a system of interacting bosons on a 
lattice in the presence of a smooth potential V(r) which 
varies on length scales much larger than the lattice spac- 
ing a. Within a local density approximation, the po- 
tential is equivalent to a shift in the local chemical po- 
tential /2(r) = ii — V(r), where fi is the global chemical 
potential determined by the total number of bosons in 
the system, N. In the situations of interest, shown in 
FigjT] the potential V(r) breaks the system into phase- 
separated domains of Mott-insulator (wherein interac- 
tions pin the number of bosons per site) and of con- 



densed bosons (which exhibit number fluctuation on each 
site). In what follows, we derive the equilibrium proper- 
ties of the domains, bulk energy costs for small deviations 
from equilibrium, dynamics of the condensed regions, the 
Josephson coupling between condensed regions mediated 
by a Mott-insulating interface, and detailed estimates for 
the spherically-symmetric situation illustrated in Figjlb. 
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FIG. 1: Top: the zero temperature phases of the Bose- 
Hubbard model; the dotted line cuts through the phases 
that could coexist for a fixed small value of zJ/U. Below: 
schematic (a) slowly-varying random potential and (b) har- 
monic confining potential, and subsequent Mott-superfluid 
domains. 

This system can be modeled by the Bose-Hubbard 
Hamiltonian, describing bosons whose tunneling between 
neighboring lattice sites has strength J and whose on-site 
repulsive interaction is U. For small J/U, each superfluid 
region is energetically near two Mott insulating phases, 
say of occupation no and no + 1. To describe this super- 
fluid region, we employ a pseudo-spin formulation of the 
Bose-Hubbard modeljEJ |9j [10] that truncates the Hilbcrt 
space to these two occupation numbers on each site. This 
formulation can be generalized to include more number 



states if necessary. But here, for simplicity and as real- 
ized in cold atom systems, we assume that J/U is suf- 
ficiently small to justify the truncation. The two-state 
Hubert space maps to a spin-1/2 basis on each site, i, 
with the identifications |no + l)j <-> | T)i an d \ n o)i I l)i> 
where the | | / J,)j are eigenstates of the spin oper- 
ator s z with eigenvalues ±1/2, and b\ — ^/n + lsf 
{pi = \friQ~+\s~), where b\ and bi denote bosonic cre- 
ation and annhilation operators, respectively, on the site 
i. The number operator, ftj — b\bi is related to the z 
component of the spin: hi = uq + 1/2 + sf. With this 
mapping, the Hamiltonian takes the form: 

H = -J(n +1) £ Of*! + aVafj+^Uno-jiiW, (1) 

(ij) • 

where (ij) denotes a summation over nearest-neighbor 
sites and /ij = /2(i"i). At the mean-field level, to which we 
confine ourselves in this Letter, the ground state config- 
uration has the pseudospins aligned with the local "mag- 
netic" field, B° = zJ{n + 1) [2/f, 2/?, cos 0J, where 
cos$i = {fii — Uno)/(zJ(no + 1)) with z the coordina- 
tion number of the lattice, the fields fi denote expecta- 
tion values of spin operators (e.g. /? = (sf)) and we 
have assumed « ij for nearest-neighbors. The equi- 
librium z-componcnt of the pseudospin has the value 
ffo = (1/2) cos 6i, the Mott phases correspond to com- 
plete polarization of the pseudospin along the z direc- 
tion, i.e. ff Q = ±1/2. Within the mean-field approxima- 
tion, we can thus identify the Mott-superfluid boundaries 
shown in Fig. [l]as occurring at the critical values of the 
external potential /x— V(v c ± ) = Un$±z J(n ±l), where ± 
refers to the boundary at the Mott no + 1 and uq phases 
respectively. 

In the condensed phase, a local order parameter can 
be defined as V = (fr f ) = Vn~ + 1/+ for < f z < 1/2 
and V = (6) = V«o + 1/" for -1/2 < f z < 0, corre- 
sponding to condensates of particles and holes, respec- 
tively. To first order in J/U and in the continuum 
limit, the equilibrium order parameter profile (as a func- 
tion of Jl = fj, — V(r)) follows from the normalization: 

fg = J\ — /q 2 /2. Ignoring the energy cost of varia- 
tions of /o from site to site (the Thomas-Fermi approx- 
imation), the order parameter is found to be: 



equilibrium in the Thomas- Fermi approximation is: 



ip(r) 



/ z 2 J 2 (n ±l) 2 -(A-n [/) 2 
4z 2 J 2 (n + 1) 



(2) 



This is of the same form as the Thomas-Fermi order 
parameter for a traditional condensate in an external 
potential V ex t and with interaction strength g: ipTF = 
\J(n — V ext )/g [3]. This allows us to identify the "effec- 
tive" confining potential for the superfluid between two 
Mott regions in the optical lattice system: {/J>—V ext ) e f / = 
{fi — n U) 2 /{zJ[n + 1]). The boson density in the con- 
densed phase is found from (h) — (no + 1 /2) + f z and in 



(ft) = (n + 1/2) + 



ft, — ngU 
2zJ(n + 1) 



(3) 



which smoothly interpolates between densities of n + 1 
at r+ and n at r~. 

For mesoscopic superfluid regions, the energy cost for 
deviations from equilibrium is non-negligible and is de- 
scribed by the bulk energy Eq. Within the Thomas- 
Fermi approximation, the Hamiltonian, Eq. 0, can be 
expressed in terms of / 2 (r): 

E B (N) = J i [zJ(n + l)(f z2 - 1) + (n U - u)f z ] dr, 

(4) 

where, assuming that variations in the density are over 
length scales greater than the lattice spacing, a contin- 
uum approximation has been made. In this approxima- 
tion, the Mott and superfluid regions are decoupled from 
one another and have separate contributions to the bulk 
energy of the system, E B (N) = E^ ott {N M ) + E S J(N S ), 
where Nm and N$ are the total number of particles in the 
Mott and superfluid phases, respectively. As seen above, 
in equilibrium, (described by the configuration /o(r)), 
the Mott-insulating particles can be thought of as pro- 
viding an effective potential that confines the superfluid 
particles. As appropriate to Josephson physics, one can 
consider a situation in which the superfluid region slightly 
shrinks or enlarges from its equilibrium configuration by 
transferring a small number of particles SN to or from 
the Mott region. In this situation, the bulk energy takes 
the form 



E f 



E^ ott (N MO ) + E s J(Nso) 



Q2 E Mott 



d 2 E s J 
ON 2 



(SN) 2 , (5) 



where the subscript '0' denotes equilibrium. The en- 
ergy scale for transfer of particles to the superfluid re- 
gion, Ec (often called the "capacitive energy" in refer- 
ence to Josephson physics in mesoscopic superconduc- 
tors) is defined by Eb = Ec(SN) 2 /2 and can be ex- 
plicitly calculated from Eq.(|5]) for a given external po- 
tential V(r). We observe that Eq.([5]) implies that the 
bulk-energy depends quadratically, rather than linearly, 
on the number of particles transfered. This differs from 
the result for two externally-trapped superfluids, where 
the linear contribution to the bulk energy only vanishes 
when the energy change of the two coupled superfluid 
regions is combined [7 . In the superfluid- Mott coexist- 
ing phases, the transfer of particles is a local one be- 
tween one superfluid region and the surrounding Mott 
phase, whose boundary is determined by the external 
potential and the ratio J/U and is in this sense self- 
organized. The coexisting system thus obeys an equilib- 
rium condition between the Mott and superfluid regions: 



2 



Q E Mott 



IdN, 



dE s J/dN s 



rendering the lowest or- 



der dependence on 6N quadratic as opposed to linear. 

Turning to the dynamics governing the Mott- 
superfluid system, in the pseudo-spin approximation the 
spins obey Heisenberg equations of motion. Further- 
more, in the mean- field approximation, the spins obey 
Bloch equations, d t ii = fj X Bj. To properly capture 
the Josephson coupling between supcrfluid regions, we go 
beyond the Thomas- Fermi approximation and allow spa- 
tial variations in the spin operators: J^j fj ~ zf + a 2 V 2 f 
in the continuum approximation. The resulting equa- 
tions of motion for the local order parameter can be 
used to derive the collective modes within each superfluid 
region [[TT], The equation of motion for the z-component 
of the spin system can be written in the form of a conti- 
nuity equation, dt(n) + V ■ J = 0, where 



(6) 



can be identified as the supercurrent density of bosons. 
To calculate the Josephson coupling between spatially- 
separated superfluid regions, we note that close to the 
Mott-superfluid interface, the Thomas- Fermi approxima- 
tion breaks down and the order parameter respects the 
equation (for f z < 0, i.e. the boundary with the no-Mott 
region) 

- J(n + l)a 2 VV 

f [Un zJ(n + 1)]^ + 2Jz\ip\ 2 ij}. (7) 



(A similar equation is respected close to the no + 1-Mott 
boundary where f z > 0.) This equation can be used 
to find the spatial decay of the order parameter beyond 
the Thomas-Fermi boundary of the Mott region. We re- 
mark that this Gross-Pitaevskii-type dynamics of the or- 
der parameter [3] breaks down well within the condensed 
phase. In particular, as seen above in Eqns.Q and ([3]), 
the Mott-superfluid system does not have a density of 
bosons directly proportional to the square of the order 
parameter. 

While the pseudo-spin description suffices at the Mott- 
superfluid boundary, and is in fact ideally suited to con- 
nect the magnitude of the order parameter at the bound- 
ary to its value in the bulk of the superfluid, it does not 
capture the physics deep in the Mott region between su- 
perfluid regions. In the no-Mott region, for instance, it is 
clear that we need to consider occupation numbers no + 1 
and no — 1 in addition to no- The relevant equations of 
motion for this case are easily calculated by employing a 
mean field perturbation theory [6], and one finds (given 
here to lowest non- vanishing order in ij)): 



in T dtip 



-K r V if> + aif>, 



1 



1 



n 



1 



n 



Un — fi 
and K r ■■ 



fi - U(n - 1) 



,(8) 



At the mean-field 



where /,•-• '/ ,. 
level, the Mott-superfluid boundary is captured by the 



relationship a = 0, which can be used to generate the 
Mott lobes of the Bose- Hubbard phase diagram shown in 
Fig(T| Furthermore, the equations of motion obtained by 
this approach, as required, coincide with Eq.Q close to 
the superfluid boundary (where terms of order \ip\ 3 can 
be ignored). 

We are now equipped to derive the Josephson cou- 
pling between two superfluid regions A and B, de- 
scribed by corresponding order parameters ipA^ A and 
tpB^ B , where ^a/b are rea l- Assuming a total or- 
der parameter of the form ipAS l ^ A + ij)B^ B , the con- 
tinuity equations of the two superfluid regions combine 
to give a continuity equation between the two regions: 
d t ((n) A - (n) B ) + V • J = 0, where J given by Eq.^, 
is found to have the Josephson form: 



(V 



J = 2Ja z VbVVm - V^V^s sin{(f> AB ), 



(9) 



where 4>ab = 4>A — 4>b is the relative phase between 
the superfluids. The Josephson energy is defined by 
dt(SNA^B) = — Ej sin(^4 — 4>b), where when parti- 
cles are transferred from the A region to the B region, 
8Na = —SNb = SNa^b- Ej can be explicitly calculated 
from Eq.([9| and the above continuity equation by inte- 
grating over an appropriate surface enclosing one of the 
superfluid regions. One finds that Ej is proportional to 
the overlap of the order parameters ipA and ips in the re- 
gion separating the two superfluids. In the two situations 
depicted in Fig. [l] this Josephson coupling a) behaves 
as a weak link bridging the two superfluid domains along 
the line of closest approach or b) has a radially sym- 
metric form connecting two concentric superfluid shells, 
and its evaluation can be reduced to a one-dimensional 
problem along the appropriate direction. In fact, the 
equilibrium configuration given by Eq. ([8| has a di- 
rect correspondence with the Ginzburg-Landau form for 
superconductors j2] and with the Gross-Pitaevksii form 
for a superfluid 3 trapped in a potential, given in this 
case by a(r). Hence, we can use standard techniques for 
calculating the Josephson coupling for a one-dimensional 
system [T^] and by employing the WKB approximation 
for the superfluid order parameters in the Mott region, 
we find 



Ej rj Aj exp 



- f VOM*' 
Jc 



(10) 



where Q(r') = z 2 Jaa(r'). The contour C can be eval- 
uated using the method of steepest descent and is the 
least-action path linking the two superfluids through the 
Mott-insulating barrier. Its end points correspond to the 
two turning points at the Mott-superfluid interface for 
A and B at which the function a vanishes. The con- 
stant A j depends on the precise forms of tpA and ipB- 
As in the case of condensates in free space, [12], Aj can 
be obtained by using a linearized potential approxima- 
tion and matching the boundary condition imposed at 



3 



the Mott-superfluid interface by Eq.Q. From Eq.(10l, 
a lower bound can be placed on the exponential depen- 
dence of the Josephson coupling by setting a to its max- 
imum value of l/(zJa 3 ) along the whole path C to ob- 
tain a value of exp(— ^£ab /a), where Iab is the path 
length. Strikingly, to first order, the Josephson coupling 
is dominated in an exponential manner only by the path 
length between superfluid regions which in turn is de- 
termined by the potential landscape. We remark that 



for the Bose-Hubbard system, Eq.(10l represents an ex- 



plicit derivation of the transport co-efficient postulated 
in Ref.[l on phenomenological grounds. 

To demonstrate the above formalism and to obtain 
estimates of the bulk and Josephson energies for an 
already-realized experimental system[S], we now consider 
N = 10 6 ultra-cold 87 Rb atoms in a three-dimensional 
optical lattice of spacing a — 0.43 /im (corresponding 
to a laser wavelength A — 2a), hopping parameter J — 
h x 120 Hz, and on-site repulsion U = hx 10 4 Hz confined 
by a harmonic trap V(r) — br 2 with b = h x 24 Hz//im 2 . 
This system has an inner Mott core with 2 atoms per 
site surrounded by a superfluid shell (SFA) , a Mott shell 
with n = 1 atom per site (1-Mott) and finally an outer 
superfluid shell (SFB); the Josephson coupling between 
the SFA and SFB shells is mediated through the 1-Mott 
shell. Eq.([8]) can be solved for the locations where a = 
to yield the boundaries of all the shells in the system. 
To calculate the capacitive energy, E c , one considers 
a transfer of a small number of particles from SFA to 
SFB, which leads to a change in the location of the re- 
gions' boundaries. Then, Eq.Q can be used to find the 
change in energy of the system. Linearizing the exter- 
nal potential in each superfluid region allows one to ob- 
tain the following expression for the bulk energy for shell 
systems where the coupling is through the n-Mott re- 
gion: E c = {ba 2 f/ 2 /{^VU)[{{2n + 1) 2 y/ f i/U - n)^ + 
((2n — l) 2 y/ fi/U — (n — l)) -1 ]. For the parameters de- 
tailed above, this leads to a bulk energy of Ec ~ h x 5 x 
10~ 3 Hz. The Josephson energy can be calculated us- 
ing Eq.(|9| after solving for the order parameter solutions 
near their respective boundaries using Eq.0. These so- 
lutions each display a characteristic decay length, cLa — 
y/ J(n + l)a 2 /qa and d B — y/ Jna 2 /q B , respectively [L^]. 
where qA/B — dV/dr\ rA/B is the slope of the external po- 
tential at the boundary of each superfluid shell (t A /b)- 
In terms of these quantities, the constant in Eq.((Tol takes 



the form Aj = (irJ A 2 / z)y/n(n + l)(r A r B a) / (d A d B ) 3/2 , 
where A « 0.397[12]. After a numerical integration, for 
the parameters detailed above we find Ej « Aje~ 2S w 
h x 2 x 10~ 8 Hz. Because the Josephson energy is expo- 
nentially dependent on the distance between the coupled 
superfluid regions, it may be possible to obtain a signif- 
icantly larger Josephson coupling in the case of a ran- 
dom (or pseudo-random) external potential where this 
distance could be more easily tuned. From the ener- 
gies found above, we can predict that the shell system 



will have Josephson oscillations which are in the strongly 
quantum regime (Ej <C i?c[Z]) and that the Josephson 
plasmon frequency, ujjp — y/EjEc ~ 10~ 4 Hz is quite 
small. This suggests that the system will be very slow 
(on the order of hours) to transfer particles between the 
two shells and that a phase difference initially present 
between the superfluids will remain for the duration of 
most current experiments, as can be ascertained via inter- 
ference experiments[13 . Similar estimates in disordered 
condensed matter systems, where Josephson physics is 
expected to play a major role, are in order. Furthermore, 
a more complete description of such co-existent Mott in- 
sulating and superfluid phases in inhomogeneous systems 
will need to incorporate several relevant factors such as 
going beyond mean field treatments, enlarging the trun- 
cated Hilbert space, and studying dissipative effects, for 
instance, due to quasiparticle excitations. 
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